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ABSTRACT 

We have examined the stacked radio and X-ray emission from UV-selected galaxies spectroscopically 
confirmed to he between redshifts 1.5 ^ z ^ 3.0 in the GOODS-North field to determine their average 
extinction and star formation rates (SFRs). The X-ray and radio data are obtained from the Chandra 
2 Msec survey and the Very Large Array, respectively. There is a good agreement between the X- 
ray, radio, and de-reddened UV estimates of the average SFR for our sample of z ~ 2 galaxies of 
~ 50 M0 yr~^, indicating that the locally-calibrated SFR relations appear to be statistically valid from 
redshifts 1.5 ^ z ^ 3.0. We find that UV-estimated SFRs (uncorrected for extinction) underestimate 
the bolometric SFRs as determined from the 2 — 10 keV X-ray luminosity by a factor of ^ 4.5 — 5.0 for 
galaxies over a large range in redshift from 1.0 <^ z < 3.5. 

Subject headings: galaxies: evolution — -galaxies: high-redshift — galaxies: starburst — dust, 
extinction — radio continuum: galaxies — X-rays: galaxies 



L INTRODUCTION 

Estimating global star formation rates (SFRs) of galax- 
ies typically requires using relations that can be quite un- 
certain as they incorporate a large number of assumptions 
in converting between specific and bolometric luminosities 
(e.g., assumed IMF, extinction, etc.; e.g., Adelberger & 
Steidel 2000). The varied efforts in the Great Observato- 
ries Origins Deep Survey (GOODS; Giavalisco et al. 2003) 
allow us to examine the same galaxies over a broad range 
of wavelengths to mitigate some of these uncertainties. X- 
ray, radio, and UV emission are all thought to directly 
result from massive stars and are consequently used as 
tracers of current star formation (e.g., Ranalli, Comastri, 
& Setti 2003; Condon 1992; Gallego et al. 1995). Here 
we use the X-ray, radio, and UV emission, each differently 
affected by extinction (or not at all), to determine SFRs 
of galaxies at z ~ 2. 

Observations of the QSO and stellar mass density, and 
morphological diversification all point to the epoch around 
z ^ 2 as an important period in cosmic history (e.g., Di 
Matteo et al. 2003; Chapman et al. 2003). Until recently, 
this epoch has been largely unexplored as lines used for 
redshift identification are shifted to the near-UV where de- 
tector sensitivity has been poor or to the near-IR, where 
spectroscopy is more difficult due to higher backgrounds. 
With the recent commissioning of the blue side of the Low 
Resolution Imaging Spectrograph (LRIS; Oke et al. 1995) 
on the Keck I telescope, we have for the first time been 
able to obtain spectra for large numbers of galaxies at 
these redshifts. Adding to the multi-wavelength efforts 
in the GOODS-North field, we have undertaken a pro- 
gram to identify photometric candidates in this field be- 
tween 1.5 ^ z ^ 3.0 and perform foUowup spectroscopy 
with LRIS-B (Steidel et al. 2004). This UV-selected sam- 
ple of galaxies forms the basis for our subsequent multi- 
wavelength analysis. 

Current sensitivity limits at X-ray and radio wave- 
lengths preclude the direct detection of normal star form- 
ing galaxies at z > 1.5. Nonetheless, we can use a "stack- 



ing" procedure to effectively add the emission from a class 
of objects in order to determine their average emission 
properties (e.g., Nandra et al. 2002; Brandt et al. 2001; 
Seibert, Heckman, & Meurer 2002). In this paper, we 
present a stacking analysis of the radio and X-ray emission 
from UV-selected star forming galaxies at redshifts 1.5 <^ 
z < 3.0 to cross-check three different techniques of esti- 
mating SFRs at high redshifts. Ho = 70 km s^^ Mpc^^, 
Um — 0.3, and f^A = 0.7 are assumed throughout. 

2. DATA 

The techniques for selecting galaxies at z ^ 2 are de- 
signed to cover the same range of UV properties and 
extinction to those used to select Lyman-break galaxies 
(LBGs) at higher redshifts (z ^ 3.0; Adelberger et al. 
2004). Here, we simply mention that we have two spec- 
troscopic samples at 1.5 ^ z ^ 2.5: a "BX" sample of 
galaxies selected on the expected UnGTZ colors of LBGs de- 
redshifted to 2.0 ^ z <^ 2.5; and a "BM" sample targeting 
z 1.5-2.0. (see Adelberger et al. 2004 and Steidel et. al. 
2004 for a complete description). We presently have 138 
redshifts ((z) ~ 2.2±0.3) and 48 redshifts ((z) ~ 1.7±0.3) 
in the GOODS-North BX and BM samples, respectively. 

The X-ray data are from the Chandra 2 Msec survey of 
the GOODS-North region (Alexander et al. 2003). We 
made use of their raw images and exposure maps in the 
Chandra soft X-ray band (0.5 — 2.0 keV). Dividing the 
raw image by the appropriate exposure map yields an im- 
age with the count rates corrected for vignetting, exposure 
time, and variations in instrumental sensitivity. The on- 
axis soft band sensitivity is ^ 2.5 x 10^^^ ergs cm^^ s^^ 
(3 ^). 

The radio data are from the Richards (2000) Very Large 
Array (VLA) survey of the Hubble Deep Field North 
(HDF-N), reaching a 3 cr sensitivity of ~ 23/xJy beam"-'^ 
at 1.4 GHz. The final naturally- weighted image has a pixel 
size of 0'.'4 and resolution of 2'.'0, with astrometric accuracy 
of < 0'.'03. 
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3. STACKING PROCEDURE 

We divided the spectroscopic data into subsets based 
on selection (BX and BM) and rcdshift, removed sources 
with matching X-ray or radio counterparts within 2'.'5 (or 
sources whose apertures are large enough to contain emis- 
sion from a nearby extended X-ray or radio source), and 
stacked galaxies in these subsets. Four of the removed 
x-ray/radio sources are detected at 850 /im with SCUBA. 

The X-ray data were stacked using the following pro- 
cedure. We added the flux within apertures randomly 
dithered by 0'.'5 at the positions of the galaxies (targets) 
in the X-ray images to produce a signal. Similarly-sized 
apertures were randomly placed within 5" of the galaxy 
positions to sample the local background near each galaxy 
while avoiding known X-ray sources. This was repeated 
1000 times to estimate the mean signal and background. 
The Chandra PSF widens for large angles from the av- 
erage pointing (off- axis angle), and we fixed the aperture 
sizes to the 50% encircled energy (EE) radii (Feigelson 
et al. 2002) for off'-axis angles > 6'. Background included 
at large off-axis angles becomes significant due to increas- 
ing aperture size and this can degrade the total stacked 
signal. Consequently, we only stacked galaxies within the 
off-axis angle that results in the highest S/N (this varies 
for each subsample, from 6 to 8'). Including all sources in 
the stack reduces the S /N but does not aflFect the absolute 
flux value. For sources < 6' from the pointing center, the 
50% EE radius falls below 2'.' 5, and we adopted a flxed 
2'.' 5 radius aperture to avoid the possibility of placing an 
aperture off a target as a result of dithering or astromet- 
ric errors which are ^ 0'.'3 for sources very close to the 
average pointing. Stacking was performed on both the 
raw and normalized images to calculate the S/N and total 
count rate, respectively. Aperture corrections were applied 
to the raw counts and count rate. The conversion between 
coimt rate and flux was determined by averaging the count 
rate to flux conversions for the 74 optically bright X-ray 
sources in the catalogs of Alexander et al. (2003; Table 7) 
which are assumed to have a photon index of F ~ 2.0, typ- 
ical of the X-ray emission from star forming regions (e.g., 
Kim, Fabbiano, & Trinchieri 1992; Nandra et al. 2002), 
and incorporate corrections for the QE degredation of the 
ACIS-I chips. In converting flux to rest-frame luminosity, 
we assume F = 2.0 and a Galactic absorption column den- 
sity of iVn = 1.6 X 10^0 cm-2 (Stark et al. 1992). Uncer- 
tainties in flux and luminosity are dominated by Poisson 
noise and not the dispersion in measured values for each 
stacking repetition, so we assume the former. 

To stack the radio data, we extracted sub-images at the 
locations of the targets from the mosaicked radio data of 
Richards (2000). These were corrected for the primary 
beam attenuation of the VLA with a maximum gain cor- 
rection of 15%, coadded using a I/ct^ weighted average to 
produce a stacked signal with maximal SjN ^ 4.5, and 
smoothed by 1'.'5. The integrated flux density, 5'i,4ghz, 
and error were computed from the standard AIPS^ task 
JMFIT using an elliptical gaussian to model the stacked 
emission. We assume a synchrotron spectral index of 
7 = —0.8, typical of the non-thermal radio emission from 

^AIPS is the Astronomical Image Processing System software 
package written and supported by the National Radio Astronomy 
Observatory. 



star forming galaxies (Condon 1992). Results of the X-ray 
and radio stacks for various subsamples are presented in 
Table 1. Four subsamples contain too few sources to yield 
a robust estimate of the stacked radio flux. 

4. RESULTS AND DISCUSSION 
4.1. SFR Estimates 

The relations established at ^; = to convert luminosi- 
ties to SFRs for our ^ ~ 2 sample are adopted from the 
following sources: Kennicutt (1998a,b) for conversion of 
the 1500 - 2800 A luminosity; Ranalli et al. (2003) for 
the 2—10 keV luminosity; and Yun, Reddy, & Condon 
(2001) for the 1.4 GHz luminosity. These relations must 
be used with caution when applied to individual galaxies 
given uncertainties in the SFR relations (e.g., burst age, 
IMF) as well as the factor of ^ 2 dispersion in the cor- 
relations between different specific luminosities. However, 
they should yield reasonable results when applied to an 
ensemble of galaxies, as we have done here. 

Table 2 shows the SFR estimates based on the 2 — 10 keV 
("SFRx"), 1.4 GHz ("SFRi.4 ghz"), and UV ("SFRuv") 
luminosities, with typical error of ~ 20%. We approxi- 
mate the UV luminosity by using the 1600 A rest-frame 
flux for all samples except the highest redshift bin sample 
(2.5 < z <^ 3.0) where we use the 1800 A rest-frame fiux. 
UV-estimated SFRs were corrected for extinction using 
the observed G — TZ colors, a spectral template assuming 
constant star formation for > 10^ yr (after which the UV 
colors are essentially constant), and applying the redden- 
ing law of Calzetti et al. (2000) and Meurer, Heckman, 
& Calzetti (1999). We created 4 additional subsamples 
of galaxies according to de-reddened UV-estimated SFR, 
also shown in Table 2. 

4.2. Stacked Galaxy Distribution and AGN 

Stacking only indicates the average emission properties 
of galaxies, not their actual distribution, and the observed 
signal may result from a few luminous sources lying just 
below the detection threshold. To investigate this, we plot 
the average distribution in counts for the sample of 147 
stacked spectroscopic galaxies (Figure 1). Much of the 
high-end tail of the distribution results from random pos- 
itive fluctuations. Only 3 sources consistently had > 7 
counts. Removing those objects whose apertures have > 6 
counts still results in a stacked signal with S/N ^ 2.5 and 
an average loss of 21 galaxies (~ 14% of the sample). It 
is therefore likely that most of the stacked galaxies con- 
tribute to the signal, particularly given their wide range 
in optical, and likely X-ray, properties. 

Contribution of low luminosity AGNs to the stacked sig- 
nal is a concern. This is a problem with most X-ray stack- 
ing analyses, but we also possess the UV spectra for our 
sources. There are two objects undetected in X-rays for 
which the UV spectra show emission lines consistent with 
an AGN. Our ability to identify AGNs from their UV spec- 
tra regardless of their X-ray properties, and having iden- 
tified only 2 such objects out of 149, suggests that sub- 
threshold luminous AGNs do not contribute significantly 
to the stacked X-ray fiux. Furthermore, UV selection bi- 
ases against the dustiest sources so we do not expect to 
find many Compton-thick AGNs in our sample. 

There are also two BM galaxies coincident with known 



Table 1 

Radio and X-ray Stacking Results 
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2.09 
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5.65 ± 0.68 


2.11 ± 0.25 


2.30 ±0.65 


5.90 ±1.66 


3.50 


BX 
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1.5 < z < 2.0 
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1.84 ±0.33 
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'^Number of galaxies stacked 
''Mean rcdshift 

■^Signal-to-noise calculated in a manner analogous to that in Nandra et al. 2002 

Average soft-band X-ray flux per object 

® Average rest-frame X-ray luminosity per object, assuming F = 2.0 and A^h = 1.6 x 10^° cm~^, for our adopted 
cosmology 

^Average integrated radio flux density per object 

^Average rest-frame 1.4 GHz luminosity per object, assuming synchrotron spectral index of 7 = —0.8 

•^Average UV luminosity computed from G and TZ magnitudes approximating the 1600 and 1800 A fluxes, respectively. 
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Star Formation Rate Estial4ter 





SFRx 


SFRr 


SER-^; 




Sample 


(M© yr-i) 


(M© yr-i) 


(M© yr-i) 


SFRx/SFR™™'- 


BX+BM 


42 


56 


50 


4.5 


BX 


42 


58 


54 


4.2 


BM 


37 




38 


4.8 


1.5 < 2 < 2.0 


37 




49 


4.3 


2.0 < 2 < 2.5 


47 




57 


4.4 


2.5 < z < 3.0 


68 




70 


4.7 


SFR^°^ < 20 Mq yr-i 


14 




11 


2.3 


20 <SFR=°^ < 40 Mq yr-i 


40 




38 


4.4 


40 <SFR=°^ < 60 M0 yr-i 


44 




48 


4.7 


SFR^°^ > 60 Mq yr-i 


72 




73 


5.1 
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radio sources that arc not detected in X-rays and arc not 
included in the stacked samples. Removing such objects 
ensures excluding radio-loud AGN that might have unas- 
suming UV and X-ray properties. For comparison, the 3 a 
radio sensitivity is sufRcient to detect SFR> 170 M0 yr^^, 
a factor of 4 higher than the median SFR of our sam- 
ple based on the X-ray or de-reddened UV SFR esti- 
mates. The stacked X-ray emission indicates a SFR of 
~ 42 M© yr~^. The on-axis soft-band flux limit implies 
a sensitivity to SFRjj 186 M© yr~^ at z ^ 2, a factor of 
4.5 higher than the average SFR for spectroscopic z ^ 2 
galaxies. Stacking the radio flux for the full spectroscopic 
sample indicates average SFRs from 33 — 70 M© yr~^ 
depending on which estimator is used: the Bell (2003), 
Condon (1992), and Yun et al. (2001) calibrations give 
low, high, and median (~ 56 M© yr~^) values, assuming 
7 = —0.8. We adopted the Yun et al. (2001) conversion 
(corrected for a binning error) as it is most relevant to the 
radio luminosity range considered here. 

4.3. Bolometric Properties of z ^ 2 Galaxies 

The UV-implied reddening indicates Ay ^ 0.5 mag and 
TVh ~ 7.5 X 10^° cm~^ assuming the Galactic calibra- 
tion (Diplas & Savage 1994). For this column density, 
absorption in the 2 — 10 keV band is negligible, and we 
therefore assiimc that SFRx is indicative of the bolometric 
SFR. In this case, we find a good agreement between the 
SFRs determined from the X-ray, radio, and de-reddcned 
UV luminosities (Lx, -t'i.4 ghz, and Luy), suggesting that 
the locally-calibrated relations between specific luminosity 
and SFR remain valid within the uncertainties at z ^ 2, 
under the caution that we cannot independently test for 
these relations as we have no direct measure of iboi- 

The (I/x) and (ii.4 ghz) of spectroscopically identified 
z 2 galaxies are comparable to those of local star- 
bursts. The X-ray/FIR relation for local galaxies (Ranalli 




Fig. 1. — Average distribution of counts for the spectroscopic 
sample. The vertical line denotes the average background count 
per aperture. The number excess at high counts (> 7) results from 
random positive fluctuations. 



et al. 2003) implies (ipiR) ^ 2.6 x lO" Lq. The stacked 
LiA GHz implies (Lfir) = 1-1 x 10" Lq (Ytm ct al. 2003). 
These estimates are similar to the FIR luminosity of lumi- 
nous infrared galaxies (LIRGs), and are expected to have 
5*850^111 ^ 0.3 mJy (e.g., Webb ct al. 2003) and would 
therefore be missing in confusion-limited SCUBA surveys 
to 2 mJy. SIRTF will have the same rest-frame 7 /im 
sensitivity to 2: ~ 2 galaxies as ISO has at 2 ~ 1 for 
LpiR ^ 5 X 10^'^ Lq galaxies (e.g., Weedman, Charman- 
daris, & Zezas 2004; Flores et al. 1999). Therefore, unhke 
the stacked averages presented here, the SIRTF data will 
be the first extinction-free tracer of the SFR distribution 
of the z ^ 2 sample as the stacked galaxies should be 
individually detected at 24 /im. 

For a fair comparison between the three redshift bins 
for 1.5 < 2; ^ 3.0, we have added back those direct X- 
ray detections in the stacks for the 1.5 < 2 < 2.0 and 
2.0 < z < 2.5 samples that would not have been detected 
if they had z > 2.5. There were no such sources with 
1.5 < z < 2.0 and only one with 2.0 < z < 2.5, increasing 
(L2-10 kev) by 2% to 2.38 x lO''^ ergs 8"^ 

The distance independent ratio SFRx/SFR^'^^'"' (Ta- 
ble 2) is similar among the selection and redshift sub- 
samples indicating that on average UV-estimated SFRs 
(uncorrected for extinction) are a factor of 4.5 times 
lower than the bolometric SFRs for galaxies between red- 
shifts 1.5 < z < 3.0. For comparison, Nandra et al. (2002) 
find this factor to be ^ 5 for both the z ^ 1 BBG and 
z ^ 3 LBG populations, and the factor is comparable to 
that of local starburst galaxies (Seibert et al. 2002). The 
attenuation computed for the BX/BM sample using the 
Calzetti et al. (2000) extinction law is similar to that 
computed from SFRx/SFR{;'^'='"-. The de-reddened UV- 
estimated SFRs (SFRyy) agree well with those predicted 
from the radio continuum for the two samples for which 
radio estimates could be obtained. Finally, we note the 
factor of ~ 5 UV attenuation is similar to that advocated 
by Steidel et al. (1999) for UV-selected samples at all 
redshifts. 

The average attenuation factor increases as the SFR in- 
creases, as shown by the last 4 subsamples in Table 2, 
and is expected if galaxies with higher SFRs have greater 
dust content on average (e.g., Adelberger & Steidel 2000). 
SFRyy follows the bolometric SFR even for low luminos- 
ity systems, indicating that the observed correlations are 
not entirely driven by only the most luminous galaxies. 

5. CONCLUSIONS 

We have made significant progress in estimating and 
comparing SFRs determined from UV, X-ray, and radio 
emission from galaxies between redshifts 1-5 <^ z 3.0, 
postulated to be the most "active" epoch for galaxy evo- 
lution. The locally-calibrated SFR relations, though un- 
certain in individual systems, appear to remain statisti- 
cally valid at high redshift. Stacking the X-ray and radio 
emission from UV-selected galaxies at 2 ~ 2 indicates that 
these galaxies have an average SFR of ^ 50 Mq yr~^ and 
an average UV attenuation factor of ^ 4.5. The prospect 
of increased radio sensitivity with the E-VLA, as well as 
X-ray campaigns in different fields to similar depth as the 
2 Msec survey in the GOODS-North field, will allow for 
a more direct probe of the radio and X-ray flux distribu- 



tion for the stacked galaxies. SIRTF/MIPS 24 data 
for the GOODS-N field will trace the dusty star formation 
in 2; ~ 2 galaxies and allow for the cross-checking of the 
results presented here. 
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